arXiv:1502.07615vl [quant-ph] 26 Feb 2015 


Generation, characterization and use of 
atom-resonant indistinguishable photon pairs 


Morgan W. Mitchell 


Abstract We describe the generation of atom-resonant indistinguishable photon 
pairs using nonlinear optical techniques, their spectral purification using atomic fil¬ 
ters, characterization using multi-photon interference, and application to quantum- 
enhanced sensing with atoms. Using either type-I or type-II cavity-enhanced spon¬ 
taneous parametric down-conversion, we generate pairs of photons in the reso¬ 
nant modes of optical cavities with linewidths comparable to the 6 MHz natural 
linewidth of the Di line of atomic rubidium. The cavities and pump lasers are 
tuned so that emission occurs in a mode or a pair of orthogonally-polarized modes 
that are resonant to the Di line, at 794.7 nm. The emission from these frequency- 
degenerate modes is separated from other cavity emission using ultra-narrow atomic 
frequency filters, either a Faraday anomalous dispersion optical filter (FADOF) with 
a 445 MHz linewidth and 57 dB of out-of-band rejection or an induced dichroism 
filter with an 80 MHz linewidth and > 35 dB out-of-band rejection. Using the type- 
I source, we demonstrate interference of photon pair amplitudes against a coherent 
state and a new method for full characterization of the temporal wave-function of 
narrow-band photon pairs. With the type-II source we demonstrate high-visibility 
super-resolving interference, a high-fidelity atom-tuned NooN state, and quantum 
enhanced sensing of atoms using indistinguishable photon pairs. 


1 Introduction 

Interference of indistinguishable photons is one of the most striking non-classical 
phenomena. When two indistinguishable photons meet at a beamsplitter, each en¬ 
tering from a different port, an interference of two-photon amplitudes leads them to 
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“coalesce” and to exit in the same direction, even though there is no force causing 
them to interact mm. This and similar multi-photon interference effects have been 
harnessed for quantum teleportation (3), entanglement swapping (4), linear optics 
quantum information processing 0, quantum-enhanced sensing m, and quantum 
simulation £3 , among other applications. A similar effect concerns the interference 
of indistinguishable photon pairs sharing a single mode against other sources of pho¬ 
ton pairs. This is the mechanism by which squeezed vacuum, which consists of pairs 
of indistinguishable photons, alters the noise distribution in an interferometer 0, an 
phenomenon central to quantum-enhanced sensing Etna and continuous-variable 
quantum information processing itTTl . 

Most of the above-mentioned applications were developed with broadband pho¬ 
tons that, due to a strong spectral mismatch, could at best interact inefficiently with 
atoms or other material systems. Generation of indistinguishable, atom-resonant 
photons is an attractive goal if we wish to interact states exhibiting non-classical in¬ 
terference with atomic quantum information processors (see chapters by Leuchs & 
Sondermann, Piro & Eschner, and Slodicka, Hetet, Hennrich & Blatt), atomic quan¬ 
tum memories (see chapter by Chuu & Du), or atomic sensors ns Karnes da. 
Using cavity-QED methods, atom-resonant indistinguishable photons can be pro¬ 
duced with exquisite control over their wave-functions (see the chapter by A. Kuhn), 
but at a high cost in system complexity. 

In this chapter we describe generation of pairs of indistinguishable photons using 
cavity-enhanced spontaneous parametric down-conversion (CESPDC) [ 0 US E) 
and extremely narrow-band optical filters flUED to select the atom-resonant com¬ 
ponent of the emission. We describe also applications: using pairs of indistinguish¬ 
able high-coherence photons, we demonstrate a full measurement of the two-photon 
temporal wave function, including both amplitude and phase (2Z|, generate an atom- 
tuned NooN state ED, and use this state to probe an atomic magnetometer, demon¬ 
strating the use of indistinguishable photon pairs for sensing of an atomic system 
beyond the standard quantum limit [23 J. 


1.1 CESPDC Sources 

Spontaneous parametric down-conversion (SPDC) is a proven method to generate 
non-classical states of light, including single-mode and two-mode squeezed states, 
entangled photon pairs, and states with multiple pairs of photons. Although photon 
production by SPDC is probabilistic, heralded single photons, heralded entangled 
states, and heralded versions of more exotic states including so-called “Schrodinger 
kitten” states, can be generated from SPDC output using single-photon detectors to 
indicate when the source has produced the desired state. Quantum states can also be 
generated “on demand” if active elements are also incorporated (see the chapter by 
Yoshikawa & Furusawa). 

A strength of SPDC is the simplicity of the parametric interaction, which can 
produce very pure correlations in frequency, polarization, and photon number. In 
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contrast, parametric interactions do not naturally produce narrowband output, and 
some engineering is required to generate SPDC output that is bandwidth-compatible 
with atomic resonance lines or other material spectroscopic features. Due to the 
strong frequency correlations, narrow-band filtering of one SPDC output can be 
used to select desired frequency content in the other output, as in the chapter by Piro 
& Eschner. Placing the SPDC process inside a resonant cavity, “cavity-enhanced 
SPDC,” (CESPDC) l24l gives a resonant enhancement to photon pairs coinciding 
with the cavity modes while suppressing those frequencies that do not fit. Cavity- 
enhancement techniques often benefit from techniques such as group-velocity mis¬ 
matching and even backward-wave phase-matching that further restrict the SPDC 
output (see the chapters by Chuu & Du, de Riedmatten & Afzelius and Zhao, Bao, 
Zhao & Pan). 

In this chapter we will describe two very similar CESPDC sources, the first using 
type-I phase matching, and thus capable of producing fully-degenerate photon pairs, 
i.e. produced in the same cavity mode and thus with the same frequency, polariza¬ 
tion, and spatial mode. The second CESPDC source uses type-II phase matching, 
and is tuned so that one //-polarized mode, and one V-polarized mode are both res¬ 
onant at half the frequency of the pump. In this way, we generate photon pairs that 
are indistinguishable in all respects except for polarization. We will sometimes re¬ 
fer to these devices as sub-threshold optical parametric oscillators (OPOs), the term 
usually used in frequency conversion and continuous-variable quantum optics for 
these same devices. 


1.2 Atomic frequency filters 


Atomic filters use atomic media, usually atomic vapors, to generate intrinsically 
narrow-band spectral features, and typically achieve transmission bandwidths from 
a few MHz to a few GHz. The “traditional” uses of atomic filters are in astronomy 
E51 . laser ranging and surveying, and daylight optical communications, where they 
are used to detect signals at specific frequencies while giving excellent blocking, 
better than 1:10 5 , over a very wide rejection band. As we will describe, these fea¬ 
tures combine well with the output of our CESPDC source, which has a linewidth 
of only a few MHz, free-spectral range (FSR) of a few hundred MHz and emission 
extending over hundreds of GHz. Proper matching of an atomic filter to a CESPDC 
source can then select a single output line while efficiently rejecting the rest of the 
SPDC output. 

While atomic media naturally present strong and narrow absorption features, 
what we require is a filter with narrowband transmission. This requires some kind 
of optical trickery, to convert the absorption resonances into transmission features. 
Equally importantly, the filter must have a strong rejection of the unwanted frequen¬ 
cies. While traditional optics, in the form of cascaded optical cavities, can in princi¬ 
ple perform this task perfectly, i.e. with arbitrarily high rejection and unit transmis¬ 
sion, this approach requires careful mode-matching and active tuning of the filters, 
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Fig. 1 Experimental setup of the type-I OPO, the FADOF filter and detection system. Symbols: 
PBS: polarizing beam splitter, AOM: acousto-optic modulator, EOM: electro-optic modulator, 
APD: avalanche photodiode, BD: calcite beam displacer, WP: Wollaston prism, TOF: time-of- 
flight analyzer, SAS: saturated absorption spectroscopy, VCO: voltage controlled oscillator, PM: 
polarization maintaining fiber, HWP: half-wave plate, QWP: quarter-wave plate, PD: photodiode 


and is ultimately limited by the quality of the cavity optics (26] [27] [28). Atomic 
filters are an interesting option for their stability, multi-mode capability, high trans¬ 
mission, and very high out-of-band rejection, limited by the quality of the polarizers. 
To date, atomic filters continue to outperform cavity filters in these figures of merit. 


2 Atom-resonant indistinguishable photon pairs in a single mode 


In this section we describe a series of experiments to generate indistinguishable pho¬ 
ton pairs in a single polarization mode, using type-I parametric downconversion. 
Subsection |2.1 1 describes the source, subsection |2.2| describes the atomic filter, sub¬ 


section |2.3| describes measurements of the spectral purity achieved, subsection f2A 
demonstrates interference of two-photon amplitudes between a narrow-band SPDC 
source and a CW laser, and subsection |2.5| applies this interference to biphoton 
wave-function measurement. Additional details can be found in (21 [ [20,22, 29). 


2.1 Type-I CESPDC source 

We use a doubly-resonant degenerate OPO li30) containing a type-I PPKTP crys¬ 
tal, phase-matched for second-harmonic generation from 794.7 nm to 397.4 nm. 
A schematic is shown in Fig. [I] A continuous wave external cavity diode laser at 
794.7 nm is stabilised at a frequency (Oq , 2.7 GHz to the red of the Rb Di line cen- 
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tre. As described in ED, an electro-optic modulator (EOM) adds sidebands to a 
saturated spectroscopy absorption signal in order to get an error signal at the desired 
frequency. We double the laser frequency by cavity-enhanced second harmonic gen¬ 
eration using a LBO crystal to generate a 397.4 nm pump beam for the OPO. The 
cavity length is stabilized to maintain a TEMOO, V-polarized cavity mode at fre¬ 
quency Q)o, thus resonating the SPDC production of indistinguishable photon pairs 
at this frequency. 

With this configuration, photon pairs are generated at all the resonance frequen¬ 
cies of the OPO cavity that fall inside the 150 GHz-wide phase matching envelope 
of the PPKTP crystal. The OPO output is thus composed of hundreds of frequency 
modes, each of 8.4 MHz bandwidth, separated by the 501 MHz FSR. We note that 
in this type-I scenario, group velocity mismatching techniques are not applicable, 
because the pairs of generated photons differ only in longitudinal mode. Filtering 
can in principle separate the desired, degenerate-mode photon pairs from the back¬ 
ground pairs, but the requirements are quite stringent: The filter must be able to 
distinguish between one cavity mode and the next, separated by the FSR. Moreover, 
to achieve a high spectral purity in the output, the filter needs a very high out-of- 
band rejection. Finally, high transmission efficiency is always desirable so as to not 
lose the photon pairs. As we shall see, these requirements are well-matched to a 
specific kind of atomic frequency filter. 


2.2 A FADOF at the Rb D\ line 

We now describe a Faraday Anomalous-Dispersion Optical Filter (FADOF) at the 
Di line of Rb at 795 nm. This line, efficiently detected with Si detectors, accessi¬ 
ble with a variety of laser technologies, and showing large hyperfine splittings, is 
a favorite for coherent and quantum optics with warm atomic vapors. Applications 
include electromagnetically-induced transparency ED, stopped light (32) , optical 
magnetometry [33, 121, laser oscillators El , polarization squeezing [35, 36 [, quan¬ 
tum memory (37) . and high-coherence heralded single photons (38, 39il. 

A FADOF consists of an atomic vapor cell between two crossed polarizers, while 
a homogeneous magnetic field along the propagation direction induces circular bire¬ 
fringence in the vapor. The crossed polarizers block transmission away from the ab¬ 
sorption line, while the absorption itself blocks resonant light. Between these, Fara¬ 
day rotation just outside the Doppler profile can give high transmission for a narrow 
range of frequencies. FADOF is simple and robust, but performance depends criti¬ 
cally on optical properties of the atomic vapor. Fortunately, first-principles modeling 
of the atomic vapor agrees very well with experiment, as shown in Fig. [2] Public- 
domain codes for calculating Rb lf2~Tl and Na (40ll4T) FADOF spectra are available. 

Similar FADOFs have been developed for several other alkali atom resonances - 
Cs D 2 El and 6S 1/2 lP 3/2 ESgl lines, Rb D 2 line (451146). Rb 5S 1/2 -► 6 P 3/2 

(47) . K (three lines) (48) . Na D lines (49), and for Ca (50) . Several of these show 
transmission above 90%, and/or linewidths below 1 GHz. Filter figures of merit are 
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Frequencyl GHz | 

Fig. 2 Features of the Rb Di FADOF spectrum for a cell of 10 cm internal length containing 
natural abundance rubidium with no buffer gas. The calculated absorption spectrum at 300K is 
shown in black, with distinct peaks due to the 85 Rb and 87 Rb transitions. The calculated spectrum 
at 365 K is shown in green, and shows very strong absorption in three regions. When this medium is 
placed between crossed polarizers in the presence of a magnetic field of B = 4.5 mT, strong Faraday 
rotation at the line-edges leads to sharp transmission peaks an points A-E. Calculated transmission 
spectra are shown in blue, experimental spectra are shown in red. The feature D, which occurs 
between two lines, is particularly advantageous, as it is both stronger (higher transmission), and 
narrower than other lines, because its tail is absorbed by the 87 Rb F = 1 F' transitions. 


shown in Table [T] An important figure of merit used in the FADOF literature is the 
Equivalent Noise Bandwidth ENBW = [max v T (v)] _1 f T(v)dv, where T (v) is the 
filter transmission versus frequency v 06). For narrowband signals in broadband 
noise, a filter achieves the signal-to-noise of an ideal filter with bandwidth ENBW. 

When quantified by the ENBW, the Rb Di line gives superior performance to 
other species and lines, due to what appears to be a fortunate accident of the hy- 
perfine splittings. For either pure 85 Rb or pure 87 Rb, the FADOF transmission at 
these field strengths shows four peaks, with the strongest ones at the extremes of the 
spectrum and with long tails. The strong 87 Rb peaks are visible as peak A and E of 
Fig. [5] The strong 85 Rb peaks include the peak D and one at —2.5 GHz, but in the 
natural-abundance vapor this latter peak is completely obscured by the 87 Rb absorp¬ 
tion. The long tail of peak D is blocked by absorption from the 87 Rb F = 1 -a F' = 1 
transition, improving the ENBW. 


2.2.1 A dual-channel FADOF 

The filter used here is a small modification of that described in Section 12.21 and in 
EEl . We take advantage of the multi-mode, imaging property of the FADOF to filter 
simultaneously two orthogonal polarizations: instead of the crossed polarizers, we 
use a beam displacer before the cell, so that the two orthogonal polarizations travel 
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Table 1 Comparison of reported FADOF transmission for different atoms and wavelengths A. 
T max : peak transmission. Bj: full-width at half-maximum bandwidth of main transmission peak. 
Bn: Equivalent-noise bandwidth (ENBW). - value not reported. Reference mo shows a K 770 nm 
FADOF curve very similar to K 766 nm. 


Atom A[nm] Ref. 



Tmax 

B t [GHz] B n [GHz] 

K 

405 

(48) 



0.93 

1.2 

6 

87 Rb 

420 

l47l 



0.98 

2.5 

5.9 

Ca 

423 

ED 



0.55 

1.5 

- 

Cs 

455 

ED 



0.96 

0.9 

3.3 

Cs 

455 

(441 



0.86 

1.5 

- 

Na 

589 

(49| 



0.85 

1.9 

5.1 

Na 

590 

(49| 



0.37 

10.5 

8.3 

K 

766 

(48| 



0.96 

0.9 

5 

K 

770 

|48l 



- 

- 

- 

Rb 

780 

(46| 



0.93 

1.3 

4.7 

Cs 

852 

142) 



0.90 

0.6 

- 

B[mT] T[K] 

Rb 

795 

(211 

18.0 

353 

0.92 

0.48 

2.1 

Rb 

795 

ED 

5.9 

378 

0.91 

1.10 

2.7 

Rb 

795 

1211 

4.5 

365 

0.71 

0.45 

1.2 

Rb 

795 

ED 

2.0 

345 

0.04 

0.32 

0.8 


along independent parallel paths in the cell. After the cell we use a Wollaston prism 
to separate the near-resonant filtered light from the unrotated one. The optical axes 
of the two polarizing elements are oriented with precision mounts, and an extinction 
ratio of 1.8 x 10 -6 is reached. 

Additionally, the setup has been supplemented with a half-waveplate placed be¬ 
fore the Wollaston prism (HWP 2 in Fig. |T]), which enables us to, in effect, turn on 
and off the filter. In the “FADOF on” condition, the waveplate axis is set parallel 
to the Wollaston axis (and thus the waveplate has no effect on the filter behaviour), 
the magnetic field is 4.5 mT and the temperature is 365 K. In the “FADOF off” 
condition, no magnetic field is applied, the temperature of the cell is also 365 K and 
HWP 2 is set to rotate the polarization by 90 degrees, in effect swapping the outputs, 
so that almost all the light is transmitted through the setup without being filtered. 

We optimized the filter using a common criterion for experiments with photon 
pairs: we maximize the ratio of coincidences due to photon pairs belonging to the 
degenerate mode to coincidences due to other photon pairs. Because of energy con¬ 
servation, the two photons in any SPDC pair will have frequencies symmetrically 
placed with respect to the degenerate mode; to prevent the pair from reaching the de¬ 
tectors, it suffices to block at least one of the photons. In terms of filter performance, 
this means that it is possible to have near-perfect filtering even with transmission in 
some spectral windows away from the degenerate mode, provided the transmission 
is asymmetrical (Fig. [3]). Using this criterion we find the optimal conditions for the 
filter performance at the field and temperature values given above. The optimum 
filter performance requires the degenerate mode that should be filtered to coincide 
with the FADOF transmission peak at a fixed frequency (2.7 GHz to the red from 






Morgan W. Mitchell 



Fig. 3 Spectral matching of FADOF to CESPDC source. Upper plot: reference transmission spec¬ 
trum of room temperature natural abundance Rb (blue), filter spectrum (black) and a mirror filter 
spectrum with respect to the degenerate cavity mode (black dashed). Red shaded regions indicate 
transmission of correlated photon pairs. Lower plot: cavity output spectrum (blue) and FADOF- 
filtered cavity spectrum (red). The degenerate cavity mode coincides with the FADOF peak. Both 
figures have the same frequency scale. 


the center of the Rb Di line). 


2.3 Spectral purification of degenerate photon pairs from type-I 
CESPDC 


We now show how such an atomic filter can separate indistinguishable, atom- 
resonant photon pairs from a much stronger broadband background of non-degenerate 
photon pairs, the natural output of a sub-threshold OPO or CESPDC source. 


Using the FADOF described in Section 1.2 we observe 70% transmission of 
the degenerate mode through the filter, simultaneous with out-of-band rejection by 
57 dB, sufficient to reduce the combined non-degenerate emission to a small frac- 
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tion of the desired, degenerate mode emission. For comparison, a recently-described 
monolithic filter cavity achieved 60% transmission and 45 dB out-of-band rejection 
l28l . We test the filter by coincidence detection of photon pairs from the OPO out¬ 
put, which provides a stringent test of the suitability for use at the single-photon 
level. We observe for the first time fully-degenerate, near atom-resonant photon 
pairs, as evidenced by correlation functions and atomic absorption measurements. 
The 96% spectral purity we observe surpasses the previous record of 94% 1391 , and 
is in agreement with theoretical predictions. 

2.3.1 Detection 

The distribution of arrival times of photons in a Hanbury-Brown-Twiss configura¬ 
tion is useful to check that the filter effectively suppresses the non-degenerate modes 
of the type-I OPO described in the previous section. We collect the OPO output in a 
polarization maintaining fiber and send it through the filter setup. The filtered light 
is then coupled into balanced fiber beam splitters that send the photons to avalanche 

that allows us 



Since we are using single photon detectors, we need to reduce as much as pos¬ 
sible the background due to stray light sources in the setup. The main source of 
background light is the counter-propagating beam that we inject in the OPO in or¬ 
der to lock the cavity length to be resonant at CDo- We solve this problem using a 
chopped lock: the experiment switches at 85 Hz between periods of data acquisi¬ 
tion and periods of stabilization. During periods of data acquisition, the AOM is off, 
and thus no locking beam is present. During periods of stabilization, the AOM is 
on, and an electronic gate circuit is used to block electronic signals from the APDs, 
preventing recording of detections due to the locking beam photons. In addition, the 
polarization of the locking beam is orthogonal to that of the OPO output. 

2.3.2 Effect of filtering on arrival time distribution 

The relative arrival-time distribution for photon pairs produced by the CESPDC 
source is a Dirac comb, with a separation given by the cavity round-trip time, times a 
double exponential with a time constant given by the cavity ring-down time BLED- 
Here the cavity round-trip time is slightly less than 2 ns, whereas the resolution of 
our time-to-digital converter is 1 ns. Evidence for the comb structure is visible in 
the measured distribution, shown in Fig. [4] This agrees well with the theoretical ex¬ 
pectation, described in the Appendix, as does the 26 ns full-width at half-maximum 
(FWHM) of the double-exponential envelope. When the filter is “on,” we expect to 
see an unmodulated double exponential of the same width. While the filter blocks 
the unwanted modes, it has little reshaping effect on the degenerate mode, which is 
much narrower than, and centred on, the peak of the filter pass-band. In agreement 
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Fig. 4 Histograms of arrival time differences for FADOF off compared to theoretical model (both 
include the background due to accidental coincidences and the artefacts resulting from 1 ns reso¬ 
lution of the counting electronics). 
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Fig. 5 Histograms of the differences of arrival times of the photon pairs for FADOF on (green) and 
FADOF on with hot cell on the path (black bars at bottom of graph, look closely!). No background 
has been subtracted. 
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with these expectations we observe a double-exponential distribution with no visible 
modulation with the same FWHM width. This is shown in Fig. [5] 


2.3.3 Spectral purity 

According to the theoretical filter spectrum from l2ll . we estimate that 98% of the 
atom-resonant photon pairs come from degenerate mode (see Fig. [3]). In order to 
test how much light outside the Rubidium Di line can pass through our FADOF, we 
split the light equally between the two different polarization paths of the filter setup 
by means of a half-wave plate put before the beam displacer (HWP 1 in Fig. [TJ. 
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A natural-abundance Rb vapor cell, with 10 Torr of N 2 buffer gas and heated un¬ 
til it is opaque for resonant light, is inserted in one of the paths after the filter. The 
collisionally-broadened absorption from this cell blocks the entire FADOF transmis¬ 
sion window, allowing us to compare the arrival time histograms with and without 
the resonant component. 

The number of photons detected after passing through the hot Rb cell is compara¬ 
ble to the detector dark counts, meaning that most of the filtered light is at the chosen 
frequency (do- We define the spectral purity Ps of the FADOF as P$ = 1 — chc/cf, 
where chc for) is th e number of photon pairs which were recorded within a coin¬ 
cidence windows of 50 ns in the path with (without) the hot cell. Considering raw 
coincidences (no background subtraction), we obtain Ps =0.98, meaning that the 
filtered signal is remarkably pure, as only the 2% of the recorded pairs are out of the 
filter spectrum. This 2% agrees with measurements of the polarization extinction 
ratio with the FADOF off, i.e., it is due to technical limitations of the polarization 
optics and could in principle be improved. Knowing that 98% of the photon pairs 
transmitted through the filter within the Rb resonance come from the degenerate 
cavity mode (due to filter spectrum), we conclude that 96% of the pairs exiting the 
filter come from the degenerate mode. 


2.4 Interference of biphoton amplitudes from distinct sources 

An SPDC process naturally generates a state of the form 


\v)oc\o)+^14.10)+..., a) 

k s ,ki 

where k s and kj index the modes of the signal and idler fields, respectively, and / 
is a complex-valued function analogous to the wave function encountered in non- 
relativistic quantum mechanics. Following this analogy, one might expect that ob¬ 
servations can only indicate / up to an unobservable global phase. A bit of reflec¬ 
tion shows, however, that even the global phase of / is observable, so long as the 
two-photon part of the wave-function exists in superposition with other parts, most 
importantly the zero-photon contribution |0). 

Consider the simplest scenario, in which signal and idler are the same mode, and 
differ at most by their times of arrival t, t r . This is in fact the case when we consider 
the filtered output of the CESPDC source described above. The SPDC state becomes 

\V) 06 |0) + jdtdt'y(t,t')a\t)a\t')\0) + ..., (2) 

where \jf now plays the role of /. The phase of y/ can be made visible by interfering 
the state against another state containing both a zero-photon and a two-photon com¬ 
ponent. A natural candidate is the continuous-wave coherent state with amplitude 
a: 
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Fig. 6 Experimental setup. AOM (EOM): acousto- (electro-) optic Modulator. PBS: polarizing 
beam splitter. QWP (HWP): quarter- (half-)wave plate. PZT: piezoelectric actuator. SMF: single 
mode fiber. PMF: polarization maintaining fiber. FBS: fiber beam splitter. 


2 

|a)°c|0) + a J dta\t) |0) + y j dt dt' a\t)a\t')\Q) +... (3) 

We can imagine combining at a polarizing beamsplitter the states | iff ), with ver¬ 
tical polarization, and |a), with horizontal polarization, to obtain a state 


|*F) = Iv'V® |a) 


H 


1 + a 

,2 


jdta^t)-^ J dtdt' t i/(t,t')ay(t)ay(t') 


+ — / dtdt'a' H (t)a' H {t') +... 


|0> v <8> |0> 


H 


(4) 

(5) 


with a two-photon component 


/ 


dt dt' 


)ay(t)ciy(t ) + a H {t)a H (t ) 


|0>y<8>|0) 


H- 


( 6 ) 


Interference of these two terms can be obtained by coincidence detection in a basis 
that is neither V nor H. For example, the rate of detection of a photon pair, at times t 
and t', both in the state | (f>) = (exp[— i<j>] \H) + \V})/y/2, is found by projection onto 
\(j))® 2 to give 
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a" 


+ — ex p [2/0] 


(7) 


which clearly depends on the phase of y/. It should be noted that observing this 
interference requires both narrow-band photon pairs, so that the detection times can 
be resolved, and a stable phase relation between a and y/. To maintain this phase 
relationship, we use the OPO also as a phase-sensitive amplifier, a well-established 
technique from continuous-variable quantum optics. 

The experimental setup is shown in Fig. [6] A continuous-wave diode laser at 
794.7 nm generates both the coherent reference beam and, after being amplified 
and doubled in frequency, a 397.4 nm pump beam for the OPO, described in fl7l . 
which generates a vertically-polarized (V) squeezed vacuum state via SPDC in a pe¬ 
riodically poled KTP crystal. The cavity length is actively stabilized with a Pound- 
Drever-Hall lock, to keep one longitudinal V mode resonant at the laser frequency. 
The locking beam is H polarized, counter-propagating, and shifted in frequency 
by an acousto-optic modulator (AOM), to match the frequency of an H-polarized 
mode. The AOM RF power is chopped and the detectors are electronically gated: 
coincidence data are acquired only when the locking light is off. With these mea¬ 
sures, the contribution of locking light to the accidental coincidences background is 
minimised. 

The V-polarized squeezed vacuum is combined with the H-polarized coher¬ 
ent reference at a polarizing beamsplitter to generate a beam with co-propagating 
squeezed and reference components. A polarization transformation, chosen so that 
|0) = (exp[— i(j>]\H) + \V))/y/2 arrives to one detector, is implemented with a 
quarter- and a half-waveplate, before coupling into a polarization-maintaining fiber. 

At the fiber output, the two polarization components are filtered with the two- 
polarization FADOF described in Section |2.3| in order to isolate the squeezed 
vacuum and block with high efficiency the hundreds of non-degenerate frequency 
modes generated by the OPO. The maximum transmission frequency of this filter 
is located at 2.7 GHz to the red of the center of the rubidium Di line, and the laser 
frequency is stabilised at this particular frequency by using an integrated electro¬ 
optic modulator to add sidebands to the laser prior to the saturated absorption spec¬ 
troscopy. 

The relative phase 0 re i between the coherent and the squeezed beam is stabilized 
by a quantum noise lock: One Stokes component is detected with a balanced po- 
larimeter, and the noise power in a 3 Hz bandwidth above 500 kHz is computed 
analogically using a multiplier circuit. This signal is fed back by a servo loop to a 
piezo-electric actuator on a mirror in the pump path, to stabilize the pump phase by 
a side-of-fringe lock. A galvanometer mirror is used to switch between the single¬ 
photon counting and stabilisation setups at a frequency of ^100 Hz. The reference 
beam power is increased during the stabilization part of the cycle, to reach the shot- 
noise-limited regime optimal for detection of the squeezing and operation of the 
noise lock. Two cascaded AOMs, whose RF power is chopped synchronously with 
the galvanometer mirror, modulate the coherent reference beam power, so that it has 
high power when the light is entering the stabilisation setup and low power when 
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Fig. 7 Time-resolved interference of the two-photon amplitude from CESPDC against the two- 
photon amplitude from a coherent state, (left) arrival-time distribution, photon pair rate as a func¬ 
tion of time difference T, for photons pairs either both diagonally-polarized or both circularly- 
polarized. A r-independent coherent state contribution interferes with a CESPDC contribution that 
has a double-exponential form centred on z = 0, with a FWHM width of 26 ns. As seen in the 
two histograms, the interference can be either constructive or destructive, (right) height of the in¬ 
terference peak versus the phase (j) appearing in the two-photon state (exp[— i<j>] \H) + | V))/y/2. As 
expected for a two-photon interference, and strikingly different from single-photon interference, 
the period of oscillation is 71 . 


the photon counting part is active. The system can maintain a fixed 0 re i over several 
hours. 

Results are shown in Fig. [7] and clearly show both constructive and destructive 
interference of the two-photon wave-function against the coherent state. 


2.5 Full reconstruction of the biphoton wave-function 

The form of Eq. 0 suggests a method to measure the two-photon wave-function, 
including not only its amplitude, but also its phase. If the term \ exp[2/0] is under 
our control, it should be possible by setting this term and measuring the resulting 
coincidence rates, to infer the complex value of This idea was made precise 

in <221 . and the above setup was used to reconstruct the biphoton wave-function 
shown in Fig. [8] 

The results are consistent with a double-exponential amplitude with 26 ns full- 
width at half-maximum (FWHM), as expected for a squeezed vacuum state from an 
OPO with the 8.1 MHz FWHM bandwidth independently-measured on our system. 
The phase of Yyv > consistent with a non-zero constant value, is reconstructed with 
a statistical uncertainty that decreases with increasing | iff^J |, reaching a$ ~ ±6 de¬ 
gree near t = 0. A constant phase is expected for an ideal OPO, while a phase 
defect could signal cavity or crystal imperfections [52, 53]. The phase offset is tun- 
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Photon arrival time difference r(ns) 


Fig. 8 Squared amplitude (above) and phase (below) of the reconstructed two-photon wave func¬ 
tion for the squeezed vacuum state. The solid line shows the predicted, double exponential ampli¬ 
tude describing an ideal squeezed vacuum state from our OPO with an independently-measured 
8.1 MHz bandwidth, with the amplitude and the offset fitted to the data. Error bars show ±1<T 
statistical uncertainty assuming Poisson statistics. 


able via the side-of-fringe lock that sets the relative phase of the squeezed vacuum 
and reference, and is another indication of interference at the two-photon level. 


3 Generation of spectrally-pure, atom-resonant NooN states 


In this section we describe a series of experiments to generate atom-tuned pho¬ 


ton pairs from Type-II SPDC. Subsection 3.1 presents a motivation in terms of 


NooN states and their interest for quantum-enhanced sensing, subsection |3.2| de¬ 
scribes the SPDC source and characterization of the generated states, subsection 


3.3 describes the atomic filter, subsection 3.4 describes measurements of the spec 


tral purity achieved, and subsection 3.5 describes the application of atom-tuned 
NooN states to quantum-enhanced sensing of magnetic fields using an atomic en¬ 
semble as a sensor. 
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3.1 NooN states 


We now describe the generation of atom-tuned photon pairs of orthogonal polariza¬ 
tion, with one photon //-polarized and one V-polarized, but indistinguishable in all 
other degrees of freedom. A state of one H state can be written 

\V) = a H a V |0) = 2 (^ a L a L~ |0) (8) 

where the circular polarization modes are defined by cll = (an + iay)/V 2, cir = 
(an — iay) / y/l. When written in terms of the photon numbers hl, and hr in the L, R 
modes respectively, and using the notation |this state is 

\ ¥ ) = ~^=(\N,0) LtR +\0,N) LtR ) (9) 

with TV = 2, an example of a “NooN” state, named for the letters that appear in the 
kets |-). 

Consider a linear polarization interferometer that operates by Faraday rotation. 
This imposes a differential phase 0 between the L and R parts of the state, trans¬ 
forming the NooN state as 

^(IVO)l,* + |O,A0l,*) ->• T(|VO)l,*W^|o,a0l,«)- (10) 

Note that the phase acts TV times on the second part of the state, because it contains 
TV photons. This implies that any signal derived from this state must vary with TV0, 
implying an TV-fold acceleration of any interferometric signal. If we consider that 
0 is an unknown phase, this accelerated interference implies an TV-fold increase in 
the Fisher information t54l . allowing estimation of (/) with uncertainty 8(f) = 1/TV, 
improving upon the standard quantum limit of 50 = 1/ \/N, the best obtainable with 
non-entangled states. 

A major motivation of this work is to test the suitability of entangled states for 
quantum enhanced sensing with atoms. In atomic media, interferometric phase shifts 
are necessarily accompanied by absorption, implying deposition of energy in the 
probed medium. Absorption also degrades any quantum advantage, as described by 
recent theory (55l[56]. To further complicate matters, in real media the phase shift 
and absorption may depend on the same unknown quantity. In a trade-off of rotation 
strength versus transparency, we employ a NooN state in a 7 MHz spectral window 
detuned four Doppler widths from the nearest 85 Rb resonance. We generate this, as 
above, with a CESPDC source and an ultra-narrow atom-based filter. 
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Fig. 9 Type-II CESPDC source: PPKTP, phase-matched nonlinear crystal; KTP, compensating 
crystal; Ml-4, cavity mirrors; PBS, polarizing beam splitter; HWP, half wave plate; QWP, quarter 
wave plate; SMF, single-mode fiber; PD, photodiode. 


3.2 Type-II CESPDC source 

As laser source we use a continuous-wave (CW) diode laser, stabilized to the Di 
transition of atomic rubidium at 795 nm and then frequency doubled to generate a 
397.5 nm pump that is passed through a mode-cleaning single-mode fiber and then 
focused into the center of a 20 mm-long periodically-poled KTiOPCU (PPKTP) 
crystal in a cavity, forming the OPO (Fig. [9]). A pump beam waist of 30 pm is 
achieved with a telescope. This beam waist was chosen to be larger than the op¬ 
timum for degenerate down-conversion according to Boyd and Kleinman E3 in 
order to reduce possible effects of thermal lensing [58] and gray-tracking [59]. The 
crystal is poled for type-II degenerate down-conversion, and produces orthogonally- 
polarized signal and idler photons. Due to crystal birefringence, these photons ex¬ 
perience temporal walk-off that would, if un-compensated, render the photons tem¬ 
porally distinguishable. A second KTP crystal of the same length and crystal cut, 
but not phase-matched and rotated about the beam direction by 90°, is added to the 
long arm of the cavity in order to introduce a second walk-off equal in magnitude 
but opposite in sign ED. 

The ring cavity is formed by two flat mirrors (Ml, M2) and two concave mir¬ 
rors (M3, M4) with a radius of curvature of 100 mm. The effective cavity length of 
610 mm corresponds to a FSR of 490 MHz. This geometry provides a beam waist 
of 42 pm for the resonant down-converted beam at the center of the crystal, which 
matches the 30 pm pump beam waist. Cavity length is controlled by a piezoelec¬ 
tric transducer on mirror Ml. The output coupler M2 has a reflectivity of 93% at 
795 nm. All other cavity mirrors are highly reflecting (R > 99.9%) at 795 nm and 
highly transmitting at 397.5 nm (R < 3%) resulting in a single-pass through the 
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nonlinear crystal for the blue pump beam. The crystal end faces are AR coated for 
397.5 nm and 795 nm. The measured cavity finesse of 70 results in a cavity linewidth 
of 7 MHz. 

While the walk-off per round trip is compensated by the KTP crystal, there is an 
uncompensated walk-off of in average half a crystal-length, because of the differ¬ 
ent positions inside the PPKTP, where a photon pair could be generated. This leads 
to a remaining temporal distinguishability at the output of the cavity that is com¬ 
pletely removed by delaying the horizontally polarized photon of each pair with a 
Michelson-geometry compensator: a polarizing beam splitter, retro-reflecting mir¬ 
rors, and quarter wave-plates set to 45° introduce an adjustable delay while pre¬ 
serving spatial mode overlap. After recombination the pairs are sent through a half 
wave plate (HWP2) that together with PBS2 determines the measurement basis. 
Both output ports of PBS2 are coupled into single-mode fibers (SMF) connected to 
single photon counting modules (Perkin Elmer SPCM-AQ4C). The pulse events are 
registered and processed by coincidence electronics (FAST ComTec P7888) with a 
resolution of 1 ns. 

The OPO cavity is actively stabilized by injecting an auxiliary beam, derived 
from the same diode laser, into the cavity via the output coupler (M2). This light is 
detected in transmission by a photodiode (PD1). Frequency modulation at 20 MHz, 
applied via the laser diode current, is used to lock to the peak of the cavity transmis¬ 
sion. To eliminate the background noise caused by this auxiliary beam and to protect 
the SPCMs, the locking and measuring intervals are alternated using a mechanical 
chopper at a frequency of about 80 Hz with a duty cycle of 24%. 

A general polarization analyzer, consisting of a quarter wave plate (QWP1) fol¬ 
lowed by a half wave plate (HWP) and a polarizing beam splitter (PBS2) is used 
to determine the measurement basis as shown in Fig. [9] To generate a NOON state 
in the HIV basis another quarter wave plate (QWP2) can be added. The two output 
ports of PBS2 are coupled to single-mode fibers and split with 50:50 fiber beam 
splitters. The four outputs are connected to a set of single photon counting modules 
(Perkin Elmer SPCM-AQ4C). Time-stamping was performed by coincidence elec¬ 
tronics with a resolution of 2 ns. By considering a time window of 150 ns, which 
is longer than the coherence time of each individual photon, we can evaluate the 
coincidences between any two of the four channels. 


3.2.1 Characterization of the NooN state 

Rotating the HWP before the detection setup allows us to demonstrate the greater 
resolution available with the NooN state, as shown in Fig. [TO] (left). First, sending 
just a single polarization and detecting the rate of single-photon arrivals, we ob¬ 
serve the expected oscillations with a period, in HWP angle, of 90°. Then, sending 
NooN states and detecting in coincidence, we observe a two-fold reduction in the os¬ 
cillation period simultaneous with high visibility, due to the two-photon coherence 
of the NooN state. The sinusoidal fit function of the coincidences shows a visibility 
of 90%. 
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(c) . T (d) 



Fig. 10 Multi-photon interference and state characterization. Left panels: resolution of one- and 
two-photons states under rotations, (a) Single-photon rotation measurement showing ordinary os¬ 
cillations with a period of 90° in the HWP angle. In this measurement only the H polarized part 
of the pair-photon state was sent to the analyzer. Vertical axis shows normalized singles rate at 
the transmitted port of PBS2. (b) High-visibility super-resolving phase measurement. Normalized 
coincidence detection between reflected and transmitted port of PBS2 for a NooN state input. The 
shorter period of the coincidence counts oscillations indicates super-resolution. Right panels: re¬ 
sults of quantum state tomography to characterize the NooN state, (c) Real and (d) imaginary part 
of the polarization density matrix of the pair-photon state transformed to a two-photon NOON state 
in the H, V basis, showing a nearly equal superposition of | HH) and \VV). 


The achieved high visibility of the state is the requirement for a high-fidelity 
NOON state. We introduce another quarter wave plate (QWP2) before the analyzing 
part of the setup to create a two-photon NOON state in the HIV basis, which can 
be written l/y/2(\H\,H2) + e 2l< ^ |Vi, V^))- Since the output state of the cavity | HV) 
is already a NOON state in the circular basis | HV) = i/\/2(\L\,Li) — \R\,R 2 ))> this 
state can be transferred into a NOON state in the HIV basis by sending it through 
an additional quarter wave plate at 45 degrees. 

We also use quantum state tomography, as in (60), to measure the polarization 
density matrix of the NooN state (H. In Fig. [TO] (right) real and imaginary parts 
of the reconstructed density matrix of a NOON state are displayed. The coherence 
of the state is partly imaginary leading to a phase of 2 (j) = 0.20 between HH and 
VV components (Fig. 10 T>)), which is however of no importance in the follow¬ 
ing. The fidelity of this state with the corresponding ideal two-photon NOON state 
1 / y/2(\Hi,H2) +e 2l( ^ |Vi, V 2 )) is 99%, making the state suitable for applications such 
as phase estimation (6) . 


3.3 Induced dichroism atomic filter 


As with the Type-I OPO described earlier, only a small portion of the SPDC output 
of the Type-II OPO is atom-resonant, and a narrow-band filter is necessary to select 
this atom-resonant portion. We use an induced dichroism atomic filter, similar in 
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Fig. 11 Principle of the circular dichroism filter, (a) non-resonant photons are blocked, either by 
the first or by the second PBS. (b) resonant photons are linearly polarized after the first PBS, may 
become circularly polarized by interaction with the circularly dichroic atomic medium, and then 
have a chance of passing the second PBS. In this way, atom-resonant light can pass the filter. In 
the ideal case of a perfect circular absorber, the transmission would still only be 25%. (c) optical 
pumping, shown here on the Di line, moves atomic population toward the states F = 2, mp = 2 and 
F = 2, mp = 1. (d) from these states, the atoms cannot absorb more o+ light on the F = 2 —>■ F' = 1 
transition, creating a strong circular dichroism at this frequency. 


HWP1 




Fig. 12 Two arrangements for induced dichroism filtering, left: pumping on the Di transition with 
geometric separation of the pump and probe beams, as in j 38 J. By combining the pump and probe 
beams using pierced mirrors, the pump beam propagates in a ring-shaped beam, coaxial with the 
probe beam, to achieve optical pumping of the atoms before they reach the probe beam, but without 
geometrical overlapping that beam, right: pumping on the D 2 transition at 780 nm with overlapped 
beams, as in (55). The wavelength difference between the pump and probe allows pump light to 
be filtered using commercial interference filters (IF), and the use of YVO 4 walkoff crystals as 
polarizers allows both polarizations to be filtered and collected. In both arrangements, the pump 
and probe are counter-propagating, giving rise to sub-Doppler resonances. 


several ways to the FADOF of Section |1.2| to separate the frequency-degenerate 
output of the CESPDC from the rest of the output. The filter, described in references 
[38, 39], has an 80 MHz FWHM passband centered on Cl>nooN and > 35 dB out- 
of-band rejection, so that only atom-tuned photons are detected. A representative 
spectrum is shown in Fig. [17] 

As shown in Fig. [T5| (right), a YVO 4 crystal separates horizontally and vertically 
polarized photons by 1 mm. The polarization modes travel parallel to each other 
through a hot rubidium cell of isotopically pure 87 Rb, optically pumped by a single- 
frequency laser resonant to the F=2^F’=3 transition of the D 2 line of 87 Rb (not 
shown). Due to Doppler shifts, the optical pumping only effects a portion of the 
thermal velocity distribution, and creates a circular dichroism with a sub-Doppler 
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Fig. 13 Transmission of the circular dichroism filter, obtained with the configuration of Fig. [15| 
(left). Blank areas around 1500 MHz indicate regions where the transmission was unmeasurably 
low, due to the strong extinction of the 85 Rb transition. 


linewidth of about 80 MHz. A second YVO 4 crystal introduces a second relative 
displacement, which can re-combine or further separate the photons, depending 
on polarization. Separated photons are collected, while re-combined photons are 
blocked. A half wave plate is used to switch between the “active” configuration, 
in which only photons that change polarization in the cell are collected, and the 
“inactive” configuration, in which photons that do not change are collected. In the 
“active” configuration, the system acts as an IFM detector for polarized atoms: a 
photon is collected only if it experiences a polarization change, i.e., if it is resonant 
with the optically pumped atoms, which absorb one circular component of the pho¬ 
ton polarization state. Neighboring modes of the degenerate mode at the rubidium 
transition are already 490 MHz detuned and therefore outside of the filter linewidth 
of 80 MHz. The out-of-band extinction ratio is > 35 dB. The filter transmission is 
optimized by adjusting the overlap between pump and single-photon mode, the ru¬ 
bidium vapor temperature and the magnitude of a small orienting applied magnetic 
field. The temperature is set to 65°C, which corresponds to an atomic density of 
5 x 10 11 cm -3 . The measured filter transmission of 10.0 % for horizontal polariza¬ 
tion and 9.5 % for vertical polarization is limited by pump power and in principle 
can reach 25% l38ll . 

To avoid contamination of the single-photon mode by scattered pump light, the 
pump enters the vapor cell at a small angle and counter-propagating to the single¬ 
photon mode. Interference filters centered on 795 nm further reject the 780 nm 
pump light with an extinction ratio of > 10 5 . The measured contribution from pump 
photons is below the detectors’ dark count rate. Each output is coupled into single¬ 
mode fiber. One is detected directly on a fiber-coupled avalanche photo diode (APD, 
Perkin Elmer SPCM-AQ4C). The other is used for subsequent experiments. Photon 
detections are recorded by a counting board (FAST ComTec P7888) for later analy¬ 


sis. 
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3.4 Spectral purity measurement 




Time delay (ns) 


Fig. 14 Induced dichroism filter spectral purity measurement, (left) schematic of the setup, (right) 
arrival-time correlations with the second vapor cell cold (green) and hot (black). When the second 
vapor cell is hot, it blocks the output to the dark count level of the detectors, implying at least 94% 
atom-resonant photons. 


The CESPDC source, shown in Fig. [9} was filtered using the geometry of Fig. 
p~5] (right), to produce in-principle spectrally pure photon pairs. The geometry of 
the filter and detection setup are shown in Fig. [14] (left). To test the spectral pu¬ 
rity, one photon was detected as a herald, and the other subjected to further spectral 
filtering using a cell of Rb vapor. The cell could be maintained at room tempera¬ 
ture, causing little absorption, or at a high temperature, in which case it efficiently 
blocked resonant photons. Fig. [14] (right) shows the results: with the cold blocking 
cell, a double-exponental arrival-time distribution is observed, as expected. When 
the blocking cell is heated, the coincidences drop to the dark-count level. We esti¬ 
mate the fraction of atom-resonant photons among the heralded single photons is at 

least 94% (321. 


3.5 Quantum-enhanced sensing of atoms using atom-tuned 
NooN states. 


We now apply the atom-tuned NooN states for sensing of Faraday rotation in a 
hot atomic ensemble. Because the Faraday rotation is a resonant phenomenon, it 
is essential to have near-resonant photons for this purpose. We use a 85 Rb atomic 
spin ensemble, similar to ensembles used for optical quantum memories m and 
quantum-enhanced atom interferometry 115, [§2j 16]. Non-destructive dispersive 
measurements on these systems, used for storage and readout of quantum informa¬ 
tion or to produce spin squeezing, are fundamentally limited by scattering-induced 
depolarization (63] 61, 14, 64]. This provides an experimentally-grounded moti¬ 
vation for the idea that the number of probe photons is a limiting resource when 
measuring this system; the number of probe photons cannot be increased without 
increasing the damage to the spin ensemble. 
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Magnetometer cell 


Fig. 15 Schematic of the setup for probing of an atomic magnetometer with atom-resonant 
NooN states. Photon pairs are generated by CESPDC as in Fig. [9] passed through a heated cell 
containing isotopically-enriched 85 Rb and immersed in a longitudinal magnetic field, filtered as 
in Fig. (right), and detected with a pair of Hanbury-Brown and Twiss setups. The three possible 
coincidence outcomes: HH, HV, and VV can then be monitored as a function of the applied B 
field. 


The setup is shown schematically in Fig. 15 Narrowband NooN states at OlNooN , 
the optical frequency of the 5 2 S ly / 2 F=2 -> SH^F^l transition of the Di line of 
87 Rb, are generated by CESPDC, as described above, and sent through the ensem¬ 
ble. The ensemble of 85 Rb atoms is contained in an anti-reflection coated vapor cell 
with internal length L = 75 mm, in a temperature-controlled oven at 70° C, together 
with a 0.5% residual 87 Rb component. An applied axial magnetic field B of up to 60 
mT produces resonantly-enhanced Faraday rotation of the optical polarization. After 
leaving the vapor cell, the photons are separated in polarization, frequency filtered, 
and detected with single-photon counters. Two counters on each polarization output 
record all possible outcomes, i.e., coincidences of HH, HV, and VV polarizations. 


3.5.1 Physics of near-resonant Faraday rotation 

The cell, with an internal path of 75 mm and containing purified 85 Rb with a small 
(0.5%) admixture of 87 Rb, no buffer gas, and no wall coatings that might preserve 
polarization, is modeled as a thermal equilibrium, Doppler-broadened vapor subject 
to Zeeman shifts in the intermediate regime. The atomic structure is calculated by 
diagonalization of the atomic Hamiltonians +H^°\ where 

/7q 1S0 ^ is the energy structure of the isotope 1S0 Rb including fine-structure contri¬ 
bution, 7/^5 = gHFS J * I is the hyperfine contribution, and = b • (g/J + g/I) 
is the Zeeman contribution. All atomic parameters are taken from references 165 k 
[ 661. The matrices are numerically diagonalized to find field-dependent en¬ 

ergy eigenstates, illustrated in Fig.[l 6 | from which the complex linear optical polar- 
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I ig. 16 Energy level diagrams relevant to Faraday rotation on the Di line of Rb. a Energy levels 
of 87 Rb relevant to generation and filtering (not to scale). The frequency of the NooN state (OnooN 
is tuned to the 5 2 S 1 / 2 F=2 transition of the Di line of 87 Rb. The optical pumping 

laser of the filter, with frequency cofp, addresses the 5 2 S!/ 2 F=2 5 2 P 3 / 2 F'=3 transition of the D 2 

line of 87 Rb. The 15 nm separation from the detection wavelength allows a high extinction using 
interference filters centered on (OnooN • b Dj energy levels of the probed ensemble versus field 
strength B , showing 85 Rb levels in blue and 87 Rb levels in red. At zero field cDnooN is 1-5 GHz 
detuned from the nearest 85 Rb transition. With increasing B, the nearest 85 Rb transition moves 
closer to resonance, increasing the Faraday rotation. The Doppler-broadened absorption begins to 
overlap C0 n oo n near B = 50 mT. 


izability is calculated, including radiative damping. The complex refractive index n± 
for a± polarizations is computed including Doppler broadening and a temperature- 
dependent atom density given by the vapor pressure times the isotope fraction, and 
the transfer function for the cell is calculated from the integral of the index along 
the beam path, including the measured drop in field strength of 15% from the center 
to the faces of the cell. Transmission spectroscopy, shown in Fig. |T7j agrees well 
with theory. 

Figs. [T6| and [T7| illustrate the optical physics of this magneto-optic rotation: The 
probe photons are red detuned from the 85 Rb transitions, and experience the same 
positive contribution to the refractive index at zero field. When the B-field is applied, 
however, the nearby 85 Rb F = 3 -A F' transitions split, with the now circularly- 
polarized transitions moving closer or farther from the probe frequency in function 
of their polarization. This provides a growing refractive index contribution for one 
circular polarization, and a decreasing contribution for the other, i.e. a circular bire¬ 
fringence giving rise to polarization rotation. Due to the proximity to resonance, the 
rotation angle increases non-linearly with B until at around B = 50 mT the 85 Rb 
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Detuning (GHz) 

Fig. 17 Spectroscopic characterization of the rubidium atomic ensemble. Circles show measured 
values, curves show predictions of first-principles model (see text), a Saturated-absorption spectra 
acquired with a natural-abundance cell at room temperature, as a frequency reference. Horizontal 
axis shows detuning from the center of the D\ spectral line, b-d Transmission spectra for the 
cell containing 85 Rb plus 0.5% 87 Rb at temperatures of 22°C, 53°C and 83°C, respectively. For 
each temperature, spectra with measured field strengths (in mT) of 0, 12, 24, 37, 49, and 58 are 
shown, in order of increasing line broadening. Grey vertical line shows (0 n oo n , the probe detuning. 
This operating point gives strong Faraday rotation with low absorption over the range 0-49 mT. 
Absorption from the small residual 87 Rb component can be seen in d. For clarity, parts a-c have 
been offset vertically by 1, 0.75, and 0.5, respectively. 


F = 3 -A F' lines begin to overlap with the probe frequency and significant absorp¬ 
tion begins. 


3.5.2 Faraday rotation signals with atom-tuned NooN states 


As seen in Fig. 18 all coincidence outcomes oscillate as a function of B , with two¬ 
fold super-resolution relative to the single-photon oscillation, visible in the singles 
counts due to a small imbalance between H and V in the input state. The interfer¬ 
ence visibilities are all > 90%, well above the 33% classical limit for HH and VV 
visibility lf67l . Also shown are predicted coincidence rates lf23lL which show good 
agreement. 

An analysis of the Fisher information from these coincidence rates confirm the 
utility of entangled states for probing atomic ensembles. The NooN state achieves 
a factor 1.30±0.05 more Fisher information per photon than the standard quantum 
limit (SQL). I.e. it gives more information per photon than can be extracted using 
any state consisting of non-entangled photons. It achieves this advantage at a field of 
B — 31 mT, near the point at which absorption begins to become important. Perhaps 
more important to eventual application, the NooN state gives also an advantage in 
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Fig. 18 High-visibility super-resolving Faraday rotation probing using optical NooN states. Top 
curve: for phase reference, singles detection rate (V polarization) versus field strength B shows one 
oscillation in the range B=0 to 50 mT. Other curves: coincidence detections HH, HV, and VV show 
two oscillations in the same range (super-resolution) and high visibility. Symbols show measured 
data (no background subtracted), with ±1 a statistical uncertainties. Curves show predicted coin¬ 
cidence rates using a first-principles atomic susceptibility calculation illustrated in Fig. [16] with a 
NooN state p found by quantum state tomography. 


Fisher information per scattered photon, the figure of merit for low-damage probing. 
As described in [23 ] the NooN state gives an advantage of 1.40± 0.06 over the SQL. 

Among other things, this exercise uncovered a previously unknown feature of 
quantum-enhanced sensing in imperfect interferometers, in this case those with 
losses. Here, as one can expect for most interferometric measurements on mate¬ 
rial systems, the loss depends on the measured quantity (here B). This dependence 
makes a positive contribution to the Fisher information, offsetting the well-known 
[68]] reduction of Fisher information due to loss of photons. 


4 Conclusions 

We have described two experimental projects to generate high-spectral purity indis¬ 
tinguishable photon pairs using cavity-enhanced SPDC and extremely narrowband 
optical filters based on resonant optical effects in atomic ensembles. Combining 
these techniques, we have demonstrated bright sources of entangled photons pairs, 
with both high two-photon coherence, for example 99% fidelity with a two-photon 
NooN state, and high spectral purity, > 94% atom-resonant heralded single photons 
for a type-II source and >98% atom resonant photon pairs for a type-I source. The 
potential for interaction with atoms is clearly shown by the generation and use of 
atom-tuned NooN states to beat the standard quantum limit in non-destructive prob¬ 
ing of an atomic magnetometer. As this work goes forward, it will be interesting to 
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see how indistinguishable photon pairs can interact with other atomic systems, for 
example atomic quantum memories and atomic quantum information processors. 


Appendix: Second-order correlation functions of filtered output 


In this section we consider the second order correlation function of the field opera¬ 
tors a out in a form: 

G ( 2 ) (r) - (al ut (t)al t (t + T)a oat (t + T)a out (t)) (11) 


for multimode (unfiltered) and single-mode (filtered) output of the OPO. 

As shown by Lu et al. [69], G^ (T) describing the output of a single-mode, far- 
below-threshold OPO has the form of double exponential decay 

Gif„ ) g ,e( r ) 0ce “ |7 ’ l(ri+ ^- ( 12 > 

where the reflectivity of the output coupler is r\ = exp[— yi t], the effective reflec¬ 
tivity resulting from intracavity losses is r 2 = exp [— 721 ] and T is the cavity round- 
trip time. An ideal narrowband filter would remove all the nondegenerate cavity- 
enhanced spontaneous down-conversion CESPDC modes, reducing the G^ 2 \t) to 
G% lQ {T). This filtering effect was demonstrated in f39l for a type-II OPO and an 
induced dichroism atomic filter. 

In lf69l it is also predicted that when the filter is off, so that the output consists of 
N cavity modes, G^> ( T ) takes the form 


r (2) 

rnulti 


P) ~ G® gle (r) 


sin 2 [(21V + l)nT /t] 
(21V +1) s\r?[nT /%] 


~ G ifn gle ( 7 ’) I S(T~n t), 

n=—°o 


(13) 

(14) 


i.e., with the same double exponential decay but modulated by a comb with a period 
equal to the cavity round-trip time T. In our case the bandwidth of the output con¬ 
tains more than 200 cavity modes, and the fraction in Eq. ( [13] ) is well approximated 
by a comb of Dirac delta functions. 

The comb period of T = 1.99 ns is comparable to the ^in = 1 ns resolution of 
our counting electronics, a digital time-of-flight counter (Fast ComTec P7888). This 
counter assigns arrival times to the signal and idler arrivals relative to an internal 
clock. We take the “window function” for the ith bin, i.e., the probability of an 
arrival at time T being assigned to that bin, to be 


f(0(T) = 


1, if T G [it £>i n , (i H - 1 )tbin\ 5 

0 , otherwise. 


(15) 
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Without loss of generality we assign the signal photon’s bin as i = 0, and we 
include an unknown relative delay Tq between signal and idler due to path length, 
electronics, cabling, and so forth. For a given signal arrival time t s , the rate of idler 
arrivals in the ith bin is f dtff^ (^)^muiti(^' — t s — Tq) (U is the idler arrival time). 
This expression must be averaged over the possible t s within bin i = 0. We also in- 
elude the “accidental” coincidence rate Gacc = *bin^i^2> where R\ ,7?2 are the singles 
detection rates at detectors 1,2, respectively. The rate at which coincidence events 
are registered with i bins of separation is then 

G muiti,det« = rW dt s f(°\t s ) (dtj f®(tj ) G® lti (tj — t s - To) + G& (16) 
0)in J J 

= t g£U(»*)2- f bm dt s f^(ts + To + nr) + G^. (17) 

n =—oo ‘'bin J 0 

We take Tq is a free parameter in fitting to the data. Note that if we write 
Tq = kt\)i n + 5 then the simultaneous events fall into kth bin and 8 G [—fbin/2, % - m /2\ 
determines where the histogram has the maximum visibility due to the beating be¬ 
tween the 1 ns sampling frequency of the detection system and the 1.99 ns comb 
period. APD time resolution is estimated to be 350 ps FWHM (manufacturer’s spec¬ 
ification), i.e. significantly less than the TOF uncertainty, and is not included here. 
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